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Abstract 

Consider a unit-resistive plate in the shape of a regular polygon with In sides, in 
which even-numbered sides are wired to electrodes and odd-numbered sides are insu- 
lated. The response matrix, or Dirichlet-to-Neumann map, allows one to compute the 
currents flowing through the electrodes when they are held at specified voltages. We 
show that the entries of the response matrix of the regular 2n-gon are given by the 
differences of cotangents of evenly spaced angles, and we describe some connections 
with the limiting distributions of certain random spanning forests. 

1 Introduction 

Let n > 2 be an integer, and let V be a regular polygon with 2n sides, centered at the 
origin in C, with the midpoint of the j th side of V located on the unit circle at e ln ^ n for 
1 < j < 2n. We imagine V to be made of a unit resistive material, and we wish to determine 
how much current will flow through V when the even-numbered sides of V are wired to 
electrodes at specified voltages while the odd-numbered sides of V are kept insulated. 

We will call the even-numbered (wired) sides of V nodes and the odd-numbered (insu- 
lated) sides free edges. Current may pass through the nodes, but no current is allowed to 
enter or exit V through the free edges. We number the n nodes so that node j corresponds 
to side 2j of V . For each j, the j th node of V is wired to an electrode and held at some 
specified voltage Vj. Each voltage configuration v = (vi, . . . ,v n ) on the nodes results in 
some current output / = (Jx, . . . , /„), where Ij is the current entering V through node j. 

This problem can be rephrased in terms of the solution to a mixed Dirichlet-Neumann 
boundary problem on the domain V: Given constants v\, . . . , v n , there is a unique contin- 
uous function V (the electric potential) on V which is harmonic on the interior, equals Vj 
on side 2j, and has zero normal derivative on the free edges (see e.g., [GM051 pg. 452 or 
pg. 445]). The current Ij entering V through node j is defined to be the flux of the electric 
field Ey = — W through side 2j into V. 

Since the current output / depends linearly on the voltage configuration v, there is an 
n x n matrix A mapping the space of voltage configurations to the space of current out- 
puts. A is called the response matrix or Dirichlet-to-Neumann map for the "resistor 
network" represented by V . The Dirichlet-to-Neumann map A is defined analogously for 
more general resistive domains or for resistor networks represented by finite graphs (see e.g. 
|CdV98] . [IM98| . [CIM98J for background). Since, for any resistor network, the k th column 
of A is the action of A on the A: th standard basis vector, Aj^ will be the current entering 
the network (V in our case) through node j when node k is held at 1 volt and all other 
nodes are held at volts. Observe that Aj fe will be positive if j = k and negative otherwise. 
Although it is not obvious from this description, the response matrix for a general electrical 
network is symmetric (see |CdV98j ). In this paper we compute the response matrix A for V 
with the boundary conditions described above: 

Theorem 1.1. The response matrix A for the regular 2n-gon with alternating wired/free 
boundary conditions has entries given by 

A COt jj^y-fc+jjj-cot [g (j-fc-i)] 
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The response matrix is closely connected to the distribution of random "groves" (a gen- 
eralization of spanning trees) in a resistor network |KWQ6j . In §[2] we give some background 
on this connection with groves, which is what initially led us to study the response matrix 
A, and we briefly compare the grove model with an analogous model based on percolation. 
We then discuss a purely algebraic approach (based on groves) that Kenyon and Wilson 
[KWQ6J § 5.2] used to compute the response matrix for the regular 2n-gon in the cases 
n = 3 and n = 4. This algebraic approach, however, is not easily adapted to general n. We 
prove Theorem 11.11 in § [3l using a combination of algebraic and analytic methods. 

We mention that it is known how to compute the response matrix by using the Schwarz- 
Christoffel formula to map the polygon to a rectangle with vertical slits that correspond to 
the free edges, as shown in Figure [TJ A general 2n-gon may be conformally mapped to a 
rectangle so that one wired side (say side j) gets mapped to the top of the rectangle, the 
adjacent free sides get mapped to the sides of the rectangle, the remaining wired sides get 
mapped to intervals of the bottom side of the rectangle, and the remaining free sides get 
mapped to vertical slits. ([KW06[ § 5.2] includes a discussion of these maps.) The current 
response Aj^ is just the ratio of the length of the image of the /c th side to the height of 
the rectangle. Without going further into details, we mention that in the example of the 
regular octagon (n = 4), this approach yields 

■ dw / — 5 — dw 



rx 5 l 

Jx 4 nf=i 



w 2 -b 2 



A i-i+2 = 7^ — 2 ~ 2 where b = 



x- 



dw / =s — ; dw 



nti(^-^) l/2 J ** nti(^-^) 1/2 

and X£ = cot ((1/2 — £)tt/8). It is not at first obvious that this should simplify to 

A j!j+ 2 = 1/2 - l/y/2. 




Figure 1: On the left is the regular 2n-gon (n = 7) with alternate sides wired and free 
(insulated). The regular 2n-gon is conformally mapped to a rectangle so that one wired 
side goes to the top and the remaining wired sides go the bottom of the rectangle, while two 
free sides get mapped to the sides of the rectangle and the remaining free sides get mapped 
to vertical slits. Entries of the response matrix are given by the lengths of the wired sides 
in the rectangle divided by the height of the rectangle. A Mdbius transformation maps the 
slit rectangle to a subset of the unit disk. Each wired side and each free side is mapped to 
an arc of a circle that passes through the top of the disk, and n — 1 of the wired sides each 
get mapped to arcs of the unit circle of length 2ir/n. 

Remark. Since the current responses for the regular polygon are differences of cotangents, 
the horizontal positions of the vertical slits and the sides of the rectangle are given by 
cotangents of evenly spaced angles. If we view the slit rectangle as being embedded in the 
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upper half plane with the rectangle's bottom edge on the real axis, then there is a Mobius 
transformation of the upper half plane to the unit disk so that images of the bases of the 
vertical slits (and the two bottom corners) are evenly spaced on the unit circle (Figure [1]). 

2 The response matrix and random groves 

We give here some background on the relation between the response matrix of a graph and 
random "groves" of that graph. A grove is a forest, i.e., an acyclic collection of edges of the 
graph, such that every constituent tree of the forest contains at least one of a special set 
of distinguished vertices, which are called nodes. The upper-left panel of Figure [2] shows a 
grove on a graph with two nodes (labeled 1 and 2), and the lower- left panel shows a grove 
on a graph with three nodes (labeled 1, 2, and 3). The first grove consists of one tree, 
which contains the two nodes 1 and 2, and the path connecting the nodes is highlighted. 
The second grove consists of two trees, one of which contains nodes 1 and 3, while the other 
contains just node 2. 




|Car921 ISnTiOll IZif95] 




.2/V3-1 |KW06j l Dub06 l 

Figure 2: Examples of grove (left) and percolation (right) configurations on regions with 
two (upper) or three (lower) nodes. Random groves and percolation are both important 
special cases of the Fortuin-Kasteleyn random cluster model from statistical physics [FK72j . 
The paths connecting the nodes are highlighted, and when the grid spacing e tends to 0, 
the probability that the nodes are connected in a given way is indicated below the panels. 
(The percolation formulas are rigorous for site percolation on the triangular lattice, and 
are conjectured to hold for the bond percolation shown here.) 
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Groves of "circular planar graphs" arise naturally in combinatorics [CS04] and statisti- 
cal physics [KW06J. A graph with distinguished nodes is "circular planar" if it embeds in 
the plane and each of the nodes lies on the outer face; the nodes are numbered in a counter- 
clockwise order. If a uniformly random grove is chosen, this defines a random (noncrossing) 
partition a on the set of nodes, and we are interested in the probability distribution of this 
random partition. 

Another natural probability distribution on node partitions arises from percolation on 
the graph (where each edge is included independently with probability p); this is illustrated 
in the right-hand panels of Figure [2] for p = 1/2. There has been significant recent work in 
studying crossing events for percolation such as those shown in Figure [21 and it is interesting 
to compare the crossing event probabilities in these two models. 

For graphs with two nodes, the partition probabilities for groves follow from Kirchhoff 's 
formula [Kir90| : 

# spanning trees Pr [grove partition is 12] 1 

# two-tree forests with ~~ Prfgrove partition is 1|2] ~ R^ 2 ~ ^ 
nodes 1 and 2 disconnected 

where R\^ is the electrical resistance between nodes 1 and 2 when the graph is viewed as 
a resistor network in which each edge has unit resistance. (The relation Ai >2 = — I/-R12 
holds for two nodes only; when there are more nodes, the response matrix entries are more 
complicated functions of the pairwise resistances.) 

In the case where the graph is a fine grid restricted to an r x 1 rectangular region, 
with the two nodes being extra vertices corresponding to the left and right edges of the 
rectangle (as illustrated in Figure [2]), then i?i 2 = r, so the probability that the grove 
partition a is 12 is simply Pr[<r = 12] = 1/(1 + r). If the fine grid is restricted to a region 
different from a rectangle, with two nodes each occupying some fraction of the boundary, 
then the region may be conformally mapped to a rectangle of some aspect ratio r, with 
the two nodes getting mapped to the left and right edges. When the grid becomes very 
fine, the resistance between the nodes in the orginal domain converges to r, so in this limit 
Pr[<7 = 12] 1/(1 +r). 

The corresponding formula for critical percolation (where each edge occurs indepen- 
dently with probability 1 /2) was deduced by Cardy |Car92] (see upper-right panel of Fig- 
ure [2]) using exact but nonrigorous methods. (Ziff adapted Cardy's formula, which was 
given for the upper half plane, to rectangular regions [Zif95| . and Smirnov gave a rigorous 
proof for the related model of site percolation [SmiOlJ.) 

More recently these boundary crossing events have been studied in regions with more 
nodes ^Dub06; KW06]; this is illustrated in the lower-right (for percolation) and lower-left 
(for groves) panels of Figure [2j In the case of groves, Kenyon and Wilson [KW06j show 
how to compute the boundary partition probabilities for any number of nodes in terms of 
the entries of the response matrix A of the graph when viewed as an electrical network. It 
is convenient to abbreviate 

. Pr [grove partition is a] 

Prur) = — — — 

Pr [grove partition is 1|2| . . . |nj 

For graphs with three nodes the analogues of Kirchhoff 's formula ([1]) are 

Pr(123) = Ai )2 Ai )3 + Ai i2 A 2>3 + A 1>3 A 2i 3, 
Pr(l|23) = -A 2 , 3 , Pr(2|13) = -Ai >3> Pr(3|12) = -A lj2 , Pr(l|2|3) = 1. 1 ' 

More generally, Kenyon and Wilson [KW06] proved that for any circular planar graph with 
any number n of nodes on the outer face, if a is any noncrossing partition on {1, . . . , n}, 
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then Pr(ci) is a polynomial in the entries of the response matrix A and can be computed 
explicitly. 

It is interesting to see what these general formulas give for some nice special cases. For 
example, Dubedat |Dub06] computed the partition probabilities for (site) percolation when 
the region is the regular hexagon with 3 nodes along alternate sides of the hexagon. (The 
lower-right panel of Figure [2] illustrates this for bond percolation.) To carry out a similar 
computation for groves, we need to calculate the response matrix for the regular hexagon 
with alternating free and wired boundary conditions on its faces. 

In the introduction we mentioned an algebraic approach for computing the response 
matrix Aj j- for a regular 2n-gon. This approach [KW06, § 5.2] is based on enumerating 
trees in a random grove; we briefly describe this tree-enumeration approach together with 
some associated open problems. 

We consider random groves on a very fine grid restricted to the regular 2n-gon, with n 
nodes "wired" to every other side of the polygon (as in Figure [2]). Using the formulas from 
KW06], we may express the polynomial 

pM _V Pr[t trees in grove] , x 
^-j' Fr[n trees in grove] 

in terms of the response matrix of the graph, which in the scaling limit approaches the 
response matrix A of the regular 2n-gon. Since (by symmetry considerations) the response 
matrix for the regular 2n-gon is circulant, we may define Aij_ju = Aj± (where indices 
are identified mod n). From the above formulas ([2]) we see that, in the limit, Pz{q) = 
1 — ?>A\q + ?>A\ q 2 . For the regular octagon it turns out that 

P 4 {q) = 1 - q{4A 1 + 2A 2 ) + g 2 (6A? + 8AiA 2 + 2A|) - g 3 (4A? + 8A?A 2 + 4AiA|). 

For each grove of a circular planar graph with n nodes there is a dual grove on the dual 
graph, which contains the duals of edges not contained in the primal grove. The number of 
trees in a grove plus the number of trees in its dual grove is n + 1. Since the dual graph of 
a fine grid restricted to the regular 2n-gon with alternate wired/free boundary conditions 
is again a fine grid restricted the the regular 2n-gon with alternate free/wired boundary 
conditions, in the limit where the grid is very fine, the probability of seeing t trees in a 
random grove equals the probability of seeing n — t + 1 trees. Hence the coefficients of the 
polynomial P n {q) form a palindrome. For n = 3 this implies that 3 A 2 = 1, or Ai = — l/Vs, 
determining the response matrix. For n = 4 there are several pairs (Ai,A 2 ) which make 
Pi{q) a palindrome, but only one in which Ai and A 2 are both negative (Ai = —1/2, 
A 2 = 1/2 - l/\/2), which is enough to determine the response matrix. For general n 
the fact that P n {q) is a palindromic polynomial generates enough constraints to limit the 
coefficients of A to a zero-dimensional algebraic variety, but it is not clear that there will 
always be a unique negative solution for the A's, nor is it clear how to obtain the solution 
for general n using this algebraic approach. 

However, in the next section we explicitly compute the response matrix for the poly- 
gon V by other means, proving Theorem ll.il Combining this result with the formulas from 
[KW06] . we may write down the first few polynomials P n (<q): 

P 2 (q) = l+q 

P 3 (q) = 1 + V3q + q 2 

P 4 (q) = 1 + (1 + V2)q + (1 + V2)q 2 + q 3 

P 5 (q) = 1 + ^h + 2Vbq + (2 + V5)q 2 + \Jh + 2\/5g 3 + q 4 

P&{q) = 1 + (2 + V3)q + (3 + 2V3)q 2 + (3 + 2V% 3 + (2 + V3)q 4 + q 5 
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The constant term is of course always 1. Referring to [KW06| . the linear term is — Y2j<k Aj, fc ' 
which simplifies to cot(-7r/(2n)). It would be interesting to better understand the polynomi- 
als P n , such as for example the approximate distribution of the number of trees for large n, 
or the asymptotics of P n (l) = 1/Pr [random grove is a single tree]. 

3 Computing A 

This section is devoted to proving Theorem 11.11 In § 13.11 we identify the eigenvectors and 
eigenvalues of A using symmetry considerations, thereby finding a diagonalization of A. 
Then in § 13.21 we compute A by performing a matrix multiplication to change from the 
eigenbasis back to the standard basis. 

3.1 Complex potential and the diagonalization of A 

The key to finding the eigenvectors and eigenvalues of A is to introduce complex electric 
potentials and currents in order to exploit the symmetry of the polygon V . A complex- 
valued potential V on V can be thought of as two separate real potentials, one from the real 
part of V and the other from the imaginary part. The electric field Ey now takes values in 
C 2 rather than R 2 , and can be thought of as carrying separate real and imaginary currents. 
The current output / will now be a complex-linear function of the voltage configuration 
v. Thus we can view A as a complex-linear transformation of C n whose restriction to R n 
yields the expected real current outputs. 

Let uj = e 2m / n , and for 1 < k < n, define the voltage configuration v k by (v k )j = ^ k ■ 

Lemma 3.1. Let 1 < k < n. Then v k = (uj k , uj 2k , . . . , u/ n_1 )' c , 1) is an eigenvector of A, 
and the corresponding eigenvalue X k is the current entering V through node n under the 
voltage configuration v k . 

Proof. Let a : C n — > C n denote the function which cyclically shifts the components of a 
vector to the left: a{v\, t>2, • • • , v n -i, v n ) = (v2, v%, . . . , v n , v\). If v is a voltage configuration 
on V, o~v is the voltage configuration obtained by replacing the voltage on node j with the 
voltage on node j + 1 (where the indices are identified mod n). Because V is rotationally 
symmetric, the resulting currents will likewise be rotated clockwise by one node. That is, 

A(av) =a(Av) for all v£ C n . 

On the other hand, for each vector = {uj k ,uj 2k , . . . ,u/ n-1 ) fc , 1), we have ov^ = uj k Vk- If 
Ik = Avt is the current output resulting from Vk, then we have 

al k = a{Av k ) = A(av k ) = A(uj k v k ) = uj k T h . 

This implies (I k )j = uj jk (I k ) n , so Av k = I k = (I k )nV k - □ 

Next we compute the eigencurrent = {I k ) n by considering the harmonic conjugate 
of the potential induced by v k . Since harmonic conjugation defines a real-linear operator - 
call it 7i - on real- valued harmonic functions on V (modulo constant functions), there is a 
unique way to extend TC to a complex-linear operator on complex-valued harmonic functions 
on V . Recall that any harmonic conjugate satisfies the Cauchy-Riemann equations, 

(HV) x = -V y and (HV) y = V x , 

so that VTiV is orthogonal to VT^ at every point. The Cauchy-Riemann equations imply 
the following well-known result (see e.g., [GamOll Section III. 6]), which we shall use to 
compute the current A&: 
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Proposition 3.2. Let £1 C C be simply connected, let V : f2 —* C be harmonic, and let 
7 : [0, 1] — > 17 6e a C 1 pa£/i wit/i 7(0) = a and 7(1) = b. IfTCV is any harmonic conjugate of 
V , then the current due to V flowing from left to right across 7 is equal to 7iV{a) — 7iV(b) 
(where 'left' and 'right' are defined relative to j's orientation). 

Now let TLqV denote the harmonic conjugate of V satisfying (TCqV)(0) = iV(0), so 
that TLq = —Id, and let V k : V — ► C be the potential function induced by the voltage 
configuration v k . 

Lemma 3.3. Let £ = e tn / n (so Q 2 = u). Then HoV k takes the value —i(( 2 i~ 1 ) k on side 
2j — 1, for 1 < j < n, and has zero normal derivative on the even-numbered sides ofV (see 
Figure^). Under the potential V k , the current X k entering V through side 2n is 2sm(nk/n). 



2k 




6k 



Figure 3: Boundary conditions of V k and T~CoV k for n = 4, illustrating Lemma 13.31 Thick 
solid lines represent wired sides of V and thin dashed lines represent free sides. 



Proof. The fact that the boundary conditions of TCoVk are wired where Vk is free and vice 
versa is immediate from the Cauchy-Riemann equations: Since VT4 is oriented parallel to 
the boundary on the free edges and orthogonal to the boundary on the nodes, the reverse is 
true for VTCVk- (Boundary issues can be dealt with by using Schwarz reflection to enlarge 
the domain, moving the boundary sides to the interior.) 

Recall that rotating the voltages in the configuration v/. clockwise by one node is the 
same as multiplying them by ui k . Thus, if z is on an even-numbered side of V, then 

V k (uz) = u k V k {z). (3) 

Now for each a G C with \a\ = 1, we define the rotation operator R a {z) = az for z G C. 
With this notation, ([3]) says that the functions Vk o and u> k Vk agree on the nodes of V. 
Since they are also both continuous on V ', harmonic on the interior, and have zero normal 
derivative on the free edges, they must be equal, so (J3j) in fact holds for all z £ V. 

Since TL is linear, ([3]) implies TioVk(ioz) = ui k TioVk(z) + C for some constant C . Setting 
z = shows that C = (1 - uj k )H V k (0) = (because H V k (0) = iV k (0) = if k / n, and 
1 — u n = 0). Thus we have 

H V k (uz) = oo k n V k (z) (4) 

for all z £ V, so TCoV k has the same type of rotational symmetry as V k . 

Let a k denote the (constant) value of TtoV k on side 1. Combining Q with the fact that 
T~CoV k has alternating free/wired boundary conditions shows that the boundary conditions 
of the harmonic functions a k V k and (TCoV k ) o R^ agree, so these functions must be equal. 
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Therefore, 

HoVk = a k V k o 

H 2 V k = a k H (V k o R^) 

-Vk = a k (H V k ) o R c -i 

-V k = a 2 k V k o R ( - 2 

-V k = a 2 k uj- k V k , 

which shows that a k = ±i£ fc (note that ([3]) was used in the last step). 

If k = n then V k = 1, which implies T~CoV k = i, so a n = i = —iC, n and X n = 0. If 
1 < k < n, we use Proposition 13.21 to compute the current \ k , up to a choice of sign: 

A fc = H Q V k {e™/ 2n ) - H Q V k {e-™l 2n ) = a k - uj~ k a k = ±i( k T *C* = +2 S m(7rk/n). 

To determine the sign, note that since V k is harmonic, we have V k {z) = E[V/%(zt)] for all 
z £ P, where zt is a standard Brownian motion started at zq = z which is reflected off 
the odd sides and absorbed at the even sides, and T is the absorption time. On the even- 
numbered sides, ReVfc(z) is maximized (with value 1) on side 2n. Thus ReV k (z) < 1 for 
z in the interior of the polygon V, so the real part of the electric field on side 2n points 
into V. Hence the real part of the current entering V on side 2n must be nonnegative, so 
it is +2 sin(7r/c/n), and a k = —iC, k . □ 



3.2 Recovering A from its diagonalization 

Putting together Lemmas 13.11 and 13.31 we have 

Theorem 3.4. The response matrix A satisfies AW = WD, where W is the matrix of 
eigenvectors given by W^ k = e 2lTl j k / n f and D is the diagonal matrix of eigenvalues with 
D ktk = X k = 2sin(7r/c/n). 

We use the following lemma to compute A from Theorem 13.41 Recall that £ = e l ' K l n . 

Lemma 3.5. If m is any integer, then 

n 

C (2m±lK = -l + icoi [TL(m± 1)] . 

l=\ 

Proof. Let^bt = (i 2m ± l ) e and let f3 = YTi=\ h t- We first show that Re/3 = -1. Since 
b n _£ = —bi, the real parts of bi and b n ^.£ add to (or the real part equals if I = n/2). 
Pairing up the terms in this way, we see ReX^=l be = and hence Re/3 = Re b n = — 1. 

Now let 9 = ^ (m± 5), and define /?' = e~ ld (3. We will use the same trick as above to 
show that is imaginary, then show that this implies Im (3 = cot 9 (see Figure HJ . Note 
that 0' = ELi <*, wh ere c e = e -% = e l£(2™±l)^-i). NoW) 

Cn+1 _, = e &(2™±l)(2»-2*+l) = e ±^ e -fe(2m±l)(2^1) = _- 

so C£ and c n+ i_£ are symmetric about the imaginary axis. It follows that f3' is imaginary 
and hence arg (3 = 9 + arg j3' = 9 ± ^ . Finally, observe that 



Im/3 = Re^ • tan(arg^) = -1 • tan (9 ± § ) = cot = cot (m ± |)] 



□ 
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Figure 4: Illustration of the proof of Lemma 13.51 with n = 6 and m = 0. The first five 
powers of ( are symmetric about the imaginary axis, so the sum (3 = Yl^=i nes on the 
line Rez = —1. The first six powers of £ are symmetric about the line argz = 9 + ir/2, 
where 6 = 7r/12, so (3 also lies on this line (equivalently, j3' = e~ lS P is imaginary). 



Proof of Theorem First we observe that the matrix W in Theorem 13.41 is invertible 
with W~ 1 = ^W*. (One can easily verify this directly, or simply notice that is 
the inverse discrete Fourier transform matrix, which is unitary.) We need to show that 
the entries of the matrix WDW~ l agree with the formula for A given in Theorem 11.11 It 
follows from Theorem EES that (WD) jtk = W jjk D kjk = ( 2 i k \ k = ( 2 i k (( k - C~ k )/i. Since 
(W- 1 )^ = = C 2jk /n, we have ' 



n i n 

A j>h = (WDW~\ k = Y,{WD) u {W- l ) t , k = ~J2 C 2j£ (C £ - C £ )C 2 ' A 



in ' 
l=\ 

cot [g (j-fc + ^j-cot [g (j-fc-j)] 
n 

where the last equality follows from Lemma 13.51 □ 
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